Electron Dynamics in Slowly Varying Antiferromagnetic Texture 
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Effective dynamics of conduction electrons in antiferromagnetic (AFM) materials with slowly varying spin 
texture is developed via non-Abelian gauge theory. Quite different from the ferromagnetic (FM) case, the spin of 
a conduction electron does not follow the background texture even in the adiabatic limit due to the accumulation 
of a SU(2) non-Abelian Berry phase. Correspondingly, it is found that the orbital dynamics becomes spin- 
dependent and is affected by two emergent gauge fields. While one of them is the non-Abelian generalization of 
what has been discovered in FM systems, the other leads to an anomalous velocity that has no FM counterpart. 
Two examples are provided to illustrate the distinctive spin dynamics of a conduction electron. 
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Interplay between current and magnetization is an essential 
issue underpinning the field of spintronics [1], which is based 
primarily on the exchange interaction that couples the spins 
of conduction electrons and localized magnetic moments. In 
ferromagnetic (FM) materials with spin texture, the exchange 
interaction forces the conduction electrons to adjust their spins 
to local magnetization in the adiabatic limit where the back- 
ground texture varies slowly over space-time. This induces a 
fictitious gauge field T^v = \rt • (d^n x d v n) to the orbital 
dynamics of conduction electrons [2-4], where n = M/\M\ 
is the direction of local magnetization and denotes space- 
time derivatives. As a consequence, the interaction between 
current and the FM texture is recast into an emergent electro- 
dynamics, and many celebrated phenomena should thereby be 
easily understood. For example, the spin motive force [5] and 
the topological Hall effect [6] are well explained by the elec- 
tric and magnetic components of the Lorentz force exerted by 
the gauge field; the back reaction of the Lorentz force in turn 
provides an intuitive interpretation to the current-induced spin 
torque exerted on the magnetic texture [3, 4, 7]. 

However, the above picture apparently fails in antiferro- 
magnetic (AFM) materials, where an itinerant spin finds no 
way to follow the orientation of the local moments that alters 
within unit cells. Nevertheless, we can define a staggered or- 
der parameter n = (Ma — Mb)/2M s which can be slowly 
varying over space-time, where Ma and Mb are the alter- 
nating local moments and M s denotes their magnitudes. A 
natural question arises is how (or one step back, whether) this 
staggered order affects the dynamics of conduction electrons 
in the adiabatic limit. Despite recent theoretical [8] and ex- 
perimental [9] progress on AFM spintronics, this problem has 
never been addressed analytically. But for a full knowledge of 
spin transport in AFM materials, an analytic effective theory 
of the conduction electrons is desired. This motivates us to fill 
the need by generalizing the gauge theory that has been set up 
in FM systems. 

In this Letter, we find that in a slowly varying AFM tex- 
ture, instead of strictly following the local staggered order, the 
spin of a conduction electron obeys a linear differential equa- 
tion (5a) which reflects an internal dynamics between degen- 
erate bands in the adiabatic limit. Correspondingly, the orbital 



dynamics becomes spin-dependent and is influenced by two 
different emergent gauge fields. While one of them is shown 
to be the non-Abelian generalization of what would be respon- 
sible for the Lorentz force in the FM case, the other leads to 
an anomalous velocity which is truly new and unique to AFM 
systems. We have achieved this goal by applying the non- 
Abelian gauge theory [10-12] on a doubly degenerate band, 
where an underlying SU(2) non-Abelian Berry phase [12, 13] 
is responsible for the internal dynamics. With two examples 
demonstrating the novel properties of the spin dynamics at 
the end, this paper provides a general framework on how an 
AFM texture affects its conduction electrons. The converse 
problem, i.e., the back reaction of conduction electrons on the 
AFM background will be left for future inquiries. 

Formalism - Consider an AFM system on a bipartite lat- 
tice with the staggered order parameter n(r , t) varying slowly 
over space-time so that the system maintains local periodicity. 
n(r,t) will be treated separately from the conduction elec- 
trons and regarded as a given field. The spin of a conduction 
electron couples to the local moments by the exchange inter- 
action J(M/M S ) • a, where a is the vector of Pauli matrices 
standing for the spin operator (we set H = 1), and M flips 
sign on neighboring A and B sublattice sites. Accordingly, the 
conduction electron is described by a nearest neighbor tight- 
binding Hamiltonian locally defined around n(r,t) which is 
purely general, 



H(n(r,t)) = 



- Jn ■ a 7(fc) 
7*(fe) Jn ■ a 



(1) 



where 7(fe) = —t e is the hopping term with S con- 
necting nearest neighboring A — B sites. In general, J can be 
negative if the exchange coupling is antiferromagnetic, but we 
assume a positive J throughout this paper. 

The local band structure can be easily solved as ±e(k) 
with e(k) = y/J 2 + |7(fc)| 2 , and in the adiabatic limit we 
neglect transitions between e and —e. Without loss of gen- 
erality we will focus on the lower band — e which is dou- 
bly degenerate with the two sub-bands labeled by A and B 
(See Fig. 1), whose wave functions are \i/j a ) = e zfe,r |u a ) and 
\ijj h ) = e ik - r \u b ). The Bloch waves \u a ) = \A(k))\ t (r, £)) 
and \ub) = \B(k))\ I (r, t)) maintain local periodicity around 
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the space-time point (r, t), where | f (f, t)) and || (r, £)) are 
eigenstates of n • a. The periodic parts |A(fe)) and \B(k)) 
exhibit opposite spatial patterns, which can be schematically 
understood in Fig. 1 that different spins try to find alternating 
sites so as to be aligned with the local moments. In the semi- 
classical point of view, an individual electron is described by 
a wave packet \W) = J dkw(k)[r] a \ip a } + rjb\ipb)], where 
J dkkw 2 (k) = k c gives the center of mass momentum, and 
(Wlrl W) = r c is the center of mass position [10]. The co- 
efficients f] a and rjb reflect the relative contributions from the 
two sub-bands and are constrained by | rj a | 2 + 1 rjb 1 2 = 1 . While 
(ip a \ipb) = due to the orthogonality of the spin eigenstates, 
(A(k)\B(k)} does not vanish. This finite overlap is of central 
importance to our theory so we hereby write it as, 



£(k) = (A(k)\B(k))-- 



| 7 (fc)| 



Ve 2 - J 2 



(2) 



It reaches maximum at the Brillouin zone (BZ) center and van- 
ishes at the BZ boundary. It is obvious from Eq. (2) that the 
larger the J, the smaller the £(fe). If J tends to infinity, the 
overlap £(fe) will vanish and the two sub-bands will be effec- 
tively decoupled. In this limit, the system becomes a simple 
combination of two independent FM sub-systems. In addition, 
£(fe) is a system parameter determined by the band structure 
alone, and it is conserved since from Eq. (2) we know that the 
energy conservation e = requires £ = 0. 

To construct effective theory on the degenerate band —e, 
it is imperative to invoke the non-Abelian formalism [10-12] 
where dynamics between the A and B sub-bands introduces 
internal degree of freedom represented by the iso-spin vector 
V = {771,7/2,773} = V^rfj, where fj = [rj a ,r] b } T . Here r is 
also a vector of Pauli matrices, but the different notation is 
adopted to avoid confusions with the spin operator a. The 
equations of motion of the electron wave packet are obtained 
through the effective Lagrangian approach [10] as, 



77. [a 
-v •[« 



rr • 



kr • 



(3a) 
(3b) 
(3c) 



where r M = (t,r c ), but fc M has no temporal component, its 
spatial components are just k c . In Eqs. (3) the • and x denote 
scalar and cross products in the iso-spin vector space. Here the 
non-Abelian Berry curvatures ft are obtained from the gauge 
potentials A on the A and B sub-bands, for instance, 



i(ui\d 



^uAn + 2*4, x A v ^ 



(4a) 
(4b) 



where z, j run between a, b. The connection in BZ van- 



ishes so n k Jl = [14], but = -SVJl = d^Al is non- 



zero due to the finite overlap. Since the two sub-bands have 
opposite spins (see Fig. 1), variation of 77 implies spin mis- 
tracking with the local order n, and it worths emphasizing 
that this is due to the non-Abelian nature of the problem rather 




A band 



B band 



FIG. 1: A schematic view of the Bloch waves in the lower band. 
Sub-band A means a local spin up electron has a larger probability 
on the A sites and a smaller probability on the B sites; sub-band B 
means the opposite case. They are degenerate in energy and their 
wave functions have a finite overlap depending on the ratio of J/e. 



than any non-adiabatic process. However, the iso-spin vector 
rj itself is not gauge invariant thus does not correspond to a 
physical observable directly. We need to relate 77 to the real 
spin of the electron defined by s = (W | a \ W) (in unit of 1 /2) 
which is fully gauge invariant. After some sophisticated cal- 
culations (see [14]), we obtain our central results: 



s = (l-e)(s-n)h, 
1 



k 



-n • (Vn x s), 
1 



r = -d k e - -(s x n) ■ h d k \n^ 



(5a) 
(5b) 

(5c) 



where h = d t n + (r • V)n, and we have omitted the sub- 
script c of r c and k c for convenience. We assert that Eqs. (5) 
are the fundamental equations of motion of a conduction elec- 
tron in an AFM texture, which are represented by the joint 
evolutions of three parameters (s, r, fc). An essential charac- 
ter that distinguishes the AFM electron dynamics from its FM 
counterpart lies in Eq. (5a), from which we know that the real 
spin sofa conduction electron does not follow the order pa- 
rameter n in the adiabatic limit. Eq. (5a) is purely geometrical 
because dt can be eliminated on both sides, therefore, a path 
of n is mapped onto a path of s. The mistrack between s and 
n is due to inter- sub-band dynamics through the accumulation 
of a SU(2) non-Abelian Berry phase V exp[—i j A!^ • rdr M ] 
along the trajectory in real space [12]. Remarkably, we note 
that the connection A 7 ^ is the same as the model proposed by 
Ref. [15] if £ is identified with their parameter f(B). Subse- 
quently, the SU(2) Berry phase can be attributed to the gauge 
flux of a 't Hooft-Polyakov monopole located at the origin of 
the parameter space spanned by n. 

Discussions - An important fact about non-Abelian theory 
is that the gauge fields Q's are not gauge invariant but the 
iso-spin scalars 77 • ft appearing in Eqs. (3) are. The gauge 
freedom originates from the phase ambiguity of the local 
spin wave functions which are obtained by acting U(r,t) = 
e -i* z( f>/2 e -i* y e/2 e -icT zX /2 on the e i g enstates of a z . Here 

6(r, t) and </>(r, t) are the spherical angles specifying the di- 
rection of n, whereas x( r ^) can be chosen arbitrarily and is 
not physical, but in deriving Eqs. (5) we have fixed the gauge 
by setting \ — 0- The gauge fields (Berry curvatures) not 
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FIG. 2: (Color online) Left panel: the iso-spin vector 77 (in blue) in 
the local frame that moves with n: 77 = r]i6 + 772 </> + 77372, where 
and <fi are spherical unit vectors. Right panel: In our particular gauge, 
77 is coplanar with n and s (in red). But no matter what gauge we 
choose, S3 = 773 and s\ + s% = ^ 2 (rjf + 77I) always hold. 



only drive the evolutions of k and r, but also determine inter- 
sub-band evolutions described by the dynamics of 77, which 
has been converted to the dynamics of s. To see the intrinsic 
relationship between s and 77 more explicitly, we now derive 
from Eq. (5a) the following expression (See [14]), 



( S -n) 2 



(s x n) 2 



(6) 



which indicates that the tip of s moves on an prolate spheroid 
with the semi-major axis being the local order parameter n, 
and the semi-minor axis has length £. On the other hand, 77 is 



constrained by 77 



\Vb\ 



1 . In our particular gauge 



marked by x = 0, 77 can be pictured as a vector in the local 
frame that moves with n(r, t) hence 77 = r\\Q + 772 </> + 77371 
(Fig. 2, left), meanwhile it is coplanar with n and s (Fig. 2, 
right). For arbitrary gauges, some manipulations show that we 
always have S3 = 773 and s\ + s\ = (?(r\\ + 772). A gauge 
transformation can only change the angles between s\^ and 
771 ? 2, while 773 is gauge invariant. 

The above properties can be further understood from two 
aspects: (i) Since opposite local spins are associated with dif- 
ferent sublattice sites, we are able to infer the spin projection 
S3 of an electron along n by measuring the probability dif- 



ference 773 



\Va\ 



I r}b\ 2 on neighboring A — B sites (vice 



versa). This explains why 773 = 53 and thus it is a physical 
variable. In contrast, 771 and 772 are not physical due to the 
gauge freedom. ( ii) The reduced density matrix for the spin 
degree of freedom after tracing out the sublattice can be writ- 
ten as p s = ^(1 + a • a), thus s = Ti[p s cr] = a. On the 
spheroid we have s 2 < 1, and it results in Trp 2 < Tip s = 1, 
which suggests that the conduction electron is effectively in a 
mixed spin state. This again can be attributed to the locking 
of S3 and 773, which means the spin and sublattice degrees of 
freedom are entangled. The entanglement provides informa- 
tion on the spin orientation from the knowledge of sublattice 
so it destroys the coherence of the spin states. 

Furthermore, we turn to the orbital dynamics which is spin- 
coupled. By substituting Eq. (5a) into (5b) we get, 



k = (1 - <f)0 • n)(E + B x r), 
E 



(7) 



^(d t 6V^ - V0ft0); B = ^(W x V0), (8) 



FM spin texture 


AFM spin texture 


s — n 


s= (1 - £ 2 )(s • n)h 


U(l) Abelian Berry Phase: 


SU(2) non- Abelian Berry Phase: 


7(r) = f r ^ M dr M 


U(T)=Vexp[-if r Al-Tdr ll ] 


Dirac monopole 


't Hooft-Polyakov monopole 


k = E+ B x r 


k = (1 — £ 2 )(s • n)(E + B x r) 


r = -d k s 


r = —due — \ (s x n) • ft du In £ 



TABLE I: Comparison of the electron dynamics in FM and AFM 
spin textures. In the FM case, spin dynamics is trivial, and along a 
closed path T the electron acquires an U(l) Berry phase which can 
be regarded as the magnetic flux of a Dirac monopole. A Lorentz 
force is resulted in the orbital motion. In the AFM case, spin dynam- 
ics is non-trivial due to the mixture of the two degenerate sub-bands 
through a SU(2) non- Abelian Berry phase, which originates from 
the gauge field generated by a 't Hooft-Polyakov monopole. Conse- 
quently, a spin-dependent Lorentz force and an anomalous velocity 
appear in the orbital dynamics. 



the E and B fields here are the same as what have been dis- 
covered in FM textures, where they are responsible for the 
spin motive force [5] and the topological Hall effect [6], re- 
spectively. However, quite different from the FM case, the 
gauge charge s ■ n in Eq. (7) is spin dependent, and the fac- 
tor £ 2 results from the non-commutative term 2A!^ x A^ v in 
Eq. (4b), they both reflect the influence of the spin dynam- 
ics on the orbital motion. The parameter £ G (0, 1) plays a 
key role here: in the £ — >> 1 limit, 1 — £ 2 vanishes thus from 
Eq. (5a) and Eq. (7) we get null results s = and k = 0. In 
the other limit where £ — >• 0, the solution of Eq. (5a) reduces 
to s = ±n upon the initial condition s(0) = ±n(0), thus 
Eq. (7) reduces to the usual Lorentz force equation, by which 
the system losses the manifest non- Abelian feature and be- 
haves as two decoupled FM sub-systems. It deserves attention 
that although the spin dynamics becomes trivial in the £ —> 
limit, a subtle difference still exists in the AFM system: the 
two sub-bands have opposite gauge charges since their spins 
are of opposite directions. Therefore, the Lorentz forces ex- 
erting on the electrons from the two sub-bands are of opposite 
signs, which may lead to spin-dependent transport. 

Besides, the real space dynamics governed by Eq. (5c) also 
exhibits spin-orbit coupling through | (s x n) h dk In £ which 
is an anomalous velocity. This term is along the same direc- 
tion as dk£, so Eq. (5c) amounts to give a modified band ve- 
locity. It worth mentioning that this term is unique to AFM 
systems and has nothing to do with the anomalous velocity 
studied in FM systems and quantum Hall systems [10]. It 
comes from the fijj^, curvature that joints the real space with 
BZ, whose importance has been overlooked before. For better 
comparison, we summarize the fundamental electron dynam- 
ics of FM and AFM textures in Tab. I. 

Examples - First consider an electron passing through an 
AFM domain wall of either the Bloch type or the Neel type. 
We assume that on the incident side s is polarized along the 
staggered order parameter n. If it were a FM domain wall 
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£, the motion of s falls into two topologically distinct classes 



FIG. 3: (Color online) Spin evolutions for three different £'s when 
n(t) is moving round a cone with constant angle from the z-axis. 
Upper panels: the tip of s is confined on the intersection of the cone's 
bottom and the spheroid described by Eq. (6). Lower panels: orbits 
of the tip plotted with s x and s y . The topology of the orbits is sep- 
arated into two classes (left and right) by the critical case (middle) 
where £;? = cos 2 0/(1 + cos 2 9) and the inner cone shrinks to zero. 
The orbits are not necessarily commensurate with n. 



then the total spin rotation on the outgoing side is a topological 
quantity Acp s = +tt where +1(— 1) is the topological number 
of the (anti-) domain wall. Here for an AFM domain wall, by 
solving Eq. (5a) it is found that the spin rotation is quantized 
by a renormalized unit: Aip 3 = ±11 where 



n = 



- arctan[£ tan ^tt] if £ < \ 
arctan[£ tan£7r] if£>|, 



(9) 



in the £ -+ limit II reduces to tt, and in the £ -+ oo limit II 
vanishes. Note II only depends on the system parameter £. 

Consider a second case where n(t) is varying round a cone 
of constant semi-angel in the lab frame (see Fig. 3), which 
can be realized in a spin wave. According to Eq. (5a) we know 
that ds z = due to dn z = 0, thus the tip of s should stay in 
the bottom plane of the cone. On the other hand, we learn 
from Eq. (6) that the tip is constrained on the spheroid that 
moves with the instantaneous n(t). Therefore, the actual orbit 
traversed by the tip is contained in the intersection of the two 
constraints, and the spin s is bounded between the n-cone 
and an inner cone whose semiangle depends on the system 
parameter £. Depicted in Fig. 3, we find that the motion of s 
in the lab frame exhibits both precession and nutation, which 
is reflected by the orbits of the tip. By tuning the parameter 



COS' 



6/(1 



cos 



separated by the critical condition £;? 
or simply \ j(k)\ = J cos 6. 
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